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Abstract 

Scalar wavelets have been used extensively in the analysis of Cosmic Microwave Background 

QQ (CMB) temperature maps. Spin needlets are a new form of (spin) wavelets which were introduced 
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. in the mathematical literature by Geller and Marinucci (2008) as a tool for the analysis of spin 
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random fields. Here we adopt the spin needlet approach for the analysis of CMB polarization 
measurements. The outcome of experiments measuring the polarization of the CMB are maps of 
the Stokes Q and U parameters which are spin 2 quantities. Here we discuss how to transform these 
spin 2 maps into spin 2 needlet coefficients and outline briefly how these coefficients can be used in 
the analysis of CMB polarization data. We review the most important properties of spin needlets, 
such as localization in pixel and harmonic space and asymptotic uncorrelation. We discuss several 
statistical applications, including the relation of angular power spectra to the needlet coefficients, 
testing for non-Gaussianity on polarization data, and reconstruction of the E and B scalar maps. 
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I. INTRODUCTION 



The WMAP satellite has provided the scientific community with the highest resolution 
full-sky data of the Cosmic Microwave Background (CMB) obtained to date (1^). These 
data have allowed precise estimates of the temperature angular power spectrum Ci up to 
the third Doppler peak and thereby high precision measurements of many cosmological 
parameters. In addition to measuring the temperature fluctuations in the CMB with high 
sensitivity, WMAP has also measured the polarization of the background radiation. But the 
polarization signal is almost an order of magnitude smaller than the temperature signal and 
is therefore measured with much lower sensitivity. The low signal-to-noise level makes it 
important to have a good understanding of systematic effects, in particular correlated noise 
properties. Little is so far known about polarized foregrounds. For low signal-to-noise data, 
small errors in the understanding of systematic effects and foregrounds could lead to errors 
in the estimates of the polarization angular power spectra. 

In the near future, the situation with respect to CMB polarization will improve signifi- 
cantly. The Planck satellite will take full-sky measurements of the polarized CMB sky (j^) 
and several other ground/balloon based experiments will follow up with high sensitivity ob- 
servations in smaller regions of the sky (the CLOVER, QUIET and QUAD experiments to 
mention a few). Both ESA and NASA are planning high sensitivity full-sky satellite borne 
experiments within the next 10-20 years. 

Polarization measurements allow for estimation of three more angular power spectra, the 
Cj^ power spectrum measuring the correlation between temperature and E mode polariza- 
tion as well as Cf^ and Cf^ which are the spectra of the E and B modes of polarization. 
The TE and EE spectra at large scales are used to estimate the reionization optical depth to 
high precision, a parameter to which the temperature power spectrum is not very sensitive. 
These spectra also give independent measurements of the other cosmological parameters 
estimated from the temperature power spectrum. This does not only serve as a consistency 
check, but also improves the statistical error bars on these parameters. Finally the much 
weaker BB polarization power spectrum may at large scales be dominated by a signal aris- 
ing from a background of gravitational waves which originated from inflation. This would 
be an important confirmation of the theory of inflation, and would as well give information 
about the detailed physics of the inflationary epoch. At smaller scales, the BB spectrum is 
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dominated by EE modes which have been converted to BB modes by gravitational lensing 
of large scale structure in the more recent universe. 

Taking into account the huge amount of polarization data which will be available in the 
next 1-2 decades, as well as the important cosmological information contained in these data, 
it is clear that efficient data analysis tools will be necessary. While a large amount of data 
analysis techniques have been developed for analyzing CMB temperature data, much less 
attention has been given to analysis techniques for polarization data, also due to the lack of 
available high sensitivity observations. 

An important tool for the analysis of temperature data has been various kinds of spher- 
ical wavelets (see j^, ^, [sl and the references therein). Indeed, in the last decade, wavelet 
methods have found applications in virtually all areas where statistical methods for CMB 
data analysis are required. Just to mention a few examples, we recall foreground sub- 



traction ([6]), point source detection 



testing for non-Gaussianity (see [9|, [lO|, 111]). 



search for anisotropies 12|, ll3l] , component separation 



;H), 



cross-correlation between CMB 



and Large Sca^ 
advocated by 



e Structure Data ([15|]), and many others. Directional wavelets have been 
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17(1 . The rationale for such a widespread interest can be explained as 



follows: CMB models are best analyzed in the frequency domain, where the behaviour at 
different multipoles can be investigated separately; on the other hand, partial sky coverage 
and other missing observations make the evaluation of exact spherical harmonic transforms 
troublesome. The combination of these two features makes the time-frequency localization 
properties of wavelets most valuable. 

More recently, a new kind of wavelets has found many fruitful applications in the CMB 



literature, the so-called spherical 



needlets. Spherical needlets were introduced in the func- 



19| . and then considered for the statistical analysis of 



tional analysis literature by 

spherical random fields by ({20, l2l|), with a view to applications to the statistical analysis 
of CMB data, including estimation of the angular power spectrum, testing for Gaussianity 
and bootstrap procedures. The related literature is already quite rich: apphcations to the 
analysis of the integrated Sachs- Wolfe effect by means of cross-correlation between CMB 



and Large Scale Structure data were given in 23]; a general introduction to their use for 
CMB data analysis is given in [23(]; a discussion on optimal weight functions is provided in 
24 1; further applications include ( 25|], 26|, |27|, l28|]). The extension to the construction of 
ISl . Il9[ | to the case of non-compactly supported weight functions is provided by [29|, l30|, ISjJ , 
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where the relationship with Spherical Mexican Hat Wavelets is also investigated, leading to 
the analysis of so-called Mexican Needlets. The stochastic properties of the corresponding 



In this paper, we show how the scalar needlets which were applied to CMB temperature 
data can be extended to polarization using spin needlets. The latter were introduced in the 
mathematical literature by [34]; they can be viewed as spin-2 wavelets designed for spin-2 
fields on the sphere having the same properties when applied to CMB polarization fields as 
scalar needlets have when applied to CMB temperature fields. In particular, we shall show 
below that spin needlets enjoy both the localization and the uncorrelation properties that 
make scalar needlets a powerful tool for the analysis of CMB temperature data. The aim 
of this paper is then to introduce the spin-2 needlets, show how they can be applied to Q 
and U maps, and discuss some preliminary ideas for possible future statistical applications 
(how these maps can be reconstructed from the spin-2 needlet coefficients, how EE and 
BB spectra may be obtained directly from these coefficients, tests of non-Gaussianity and 
others) . 

The plan of the paper is as follows: in Section 2 we provide a quick review of standard 
(scalar) needlets and their main properties; in Section 3 we report the main results of [34], 
where spin needlets were first introduced and investigated from the mathematical point of 
view. In Section 4 we discuss a number of possible future applications to polarization data 
analysis; in Section 5 we provide a comparison with alternative approaches for the wavelet 
analysis of polarization data, while in Section 6 we present some preliminary evidence on 
the reconstruction properties of spin needlets. Some background material on spin spherical 
harmonics is given in the Appendix. 

II. A REVIEW OF SCALAR NEEDLETS 

To ease comparisons with the spin needlets which we will describe later, we recall very 
briefly the construction of a needlet basis. The spherical needlet (function) is defined as 



here, 7 is a direction {6, 0) on the sphere, and j is the frequency (multipole range) of the 
needlet. We use {^jk} to denote a set of cuhature points on the sphere, corresponding to 



Mexican needlet coefficients are established in [32, 





(1) 
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a given frequency j. In practice, we will identify these points with the pixel centres in the 
HEALPix pixelization scheme [ssl. The cubature weights Xjk are inversely proportional to 
the number of pixels we are actually considering (see 2^, 2l| for more details); we can take 
for simplicity Xj^. = An/Nj, where here and throughout the paper we shall take Nj to be the 
number of cubature points {C,jk} ■ Needlets can then be viewed as a combination, analogous 
to convolution, of the projection operators Ylm=-e^^rn{'y)yim{^jk) with a suitably chosen 
window function b{x) with x = i/BK The number B appearing in the argument of the func- 
tion b{x) is a parameter which defines the needlet basis, as will be discussed below. Special 
properties of b{x) ensure that the needlets enioy quasi-exponential localization properties in 



pixel space. Formally, we must ensure that 
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19|): 



Al The function is positive in the range x = B], zero otherwise; hence 
is positive in £ G [B^-\ B^+^] 

A2 the function b{x) is infinitely differentiable in (0, oo). 



A3 we have 



i=i 



1 for all i> B. 



(2) 



Condition Al can be generalized to cover functions such as xexp(— a;), thus leading to so- 



called Mexican needlets, see 



29 



30 



31| where advantages and disadvantages of such choice 



are also discussed. In the present formulation, Al ensures the needlets have bounded support 
in the harmonic domain; A2 is needed for the derivation of the localization properties, which 
we shall illustrate in the following section. Finally, A3 {the partition of unity property) is 
needed to establish the reconstruction formula ([6]). Examples of constructions satisfying 



A1-A3 are given in [22, 12^. 



We shall now recall briefly some of the general features of needlets, which are not in 
general granted by other spherical wavelet constructions. We refer for instance to 23|] for 
more details and a general introduction for a CMB readership. Briefly, we recall the following 
features: 

a) needlets do not rely on any tangent plane approximation (compare and take 
advantage of the manifold structure of the sphere; 

b) being defined in harmonic space, they are computationally very convenient, and in- 
herently adapted to standard packages such as HEALPix; 



c) they allow for a simple reconstruction formula (see ([6])), where the same needlet func- 
tions appear both in the direct and the inverse transform. This property is the same as for 
spherical harmonics but it is not shared by other wavelets systems; 

d) they are quasi-exponentially (i.e. faster than any polynomial) concentrated in pixel 
space, see ([7]) below; 

e) they are exactly localized on a finite number of multipoles; the width of this support 
is explicitly known and can be specified as an input parameter (see ([1])); 

f ) random needlet coefficients can be shown to be asymptotically uncorrelated (and hence, 
in the Gaussian case, independent) at any fixed angular distance, when the frequency in- 
creases (see I20I). This capital property can be exploited in several statistical procedures, 
as it allows one to treat needlet coefficients as a sample of independent and identically 
distributed coefficients on small scales, at least under the Gaussianity assumption. 

More precisely, random needlet coefficients are given by 

f3jk = [ W)M'f)dl (3) 
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e m=-e ^-^^^ ^ 

i ^ 

= \^k^b{—) aemYe^i^jk)- (4) 

e m=-i 

Here j denotes the frequency of the coefficient and k refers to the direction {6, (p) on the 
sky. The index k can in practice be the pixel number on the HEALPix grid. It is very 
important to stress that, although the needlets do not make up an orthonormal basis for 
square integrable functions on the sphere, they do represent a tight frame. In general, a tight 
frame on the sphere is a countable set of functions which preserves the norm; frames do not 
in general make up a basis, as they admit redundant elements. They can be viewed as the 



closest system to a basis, for a given redundancy, see [20|, l3l|, l35| for further definitions and 



discussion. In our framework, the norm-preserving property becomes 

„ 00 

jk -^^^ £=1 



where 



-E 
+ 1 ^ 



|2 



tm 



7 



is the raw angular power spectrum of the map ^(7). ([5]) suggests immediately some proce- 
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28|) and is related to a much 

(6) 



dures for angular power spectrum estimations and testing 
more fundamental result, i.e. the reconstruction formula 

which in turn is a non-trivial consequence of the careful construction leading to ([2]). Again, 
we stress that the simple reconstruction formula of ([6]) is typical of tight frames but does not 
hold in general for other wavelet systems. It is easy to envisage many possible applications 
of ([6]) when handling masked data. 



A. Localization and uncorrelation properties 

The following quasi-exponential localization property of needlets is due to Q and 
motivates their name: 

For any M = 1, 2, ... there exists a positive constant Cm such that for any point 7 G S"^ 
we have 

W'^-'^'l ^ (l + B^arccosm-O.!))" ' 
We recall that arccos(|7 — C,jk\) is just the geodesic distance on the unit sphere between 

the position 7 and the position ^^7^; (we recall in practice C,jk can be the pixel center of a 

HEALPix pixel k). ([7]) is then stating that, for any fixed nonzero geodesic distance, the value 

of ipjkil) goes to zero faster than any polynomial {quasi- exponentially) in the parameter B. 

Thus needlets achieve excellent localization properties in both the real and the harmonic 

n 

domain. In [29], (JTj) is extended to the case of a non-compactly supported but smooth b{x), 
thus covering also the Mexican needlet case (where b{x) ~ xexp(— x)). 

From the stochastic point of view, the crucial uncorrelation properties for random spher- 



ical needlet coefficients were given in 20]. More precisely, two forms of uncorrelation were 
established 

• PI Whenever \ji — ^2! > 2, we have that {PjikiPj2k2) = 0, () denoting the expected 
value 

• P2 For ji = j2, for any M > there exist a constant Cm such that 

\{(^jkf3jk')\ ^ Cm_ , s 

- (1 + arccos(|e,fc - ^,^'1))*^ ' ^ ^ 
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Property PI is a straightforward consequence of the locahzation properties for b{x) in the 
harmonic domain. Property P2 is much more surprising, and does not follow by any means 
from the localization in pixel space ([7]); indeed it is simple to provide examples of wavelet 



systems that satisfy ([7]), and still do not enjoy ([S]) (see [32j, 33|]). Both these properties hold 
for any isotropic random field, without any assumption on its distribution (i.e. Gaussianity) . 

Properties P1,P2 suggest that at high frequency, needlet coefficients can be approximated 
as a sample of identically distributed and uncorrelated (independent, in the Gaussian case) 
coefficients, and this property opens the way to a huge toolbox of statistical procedures 
for CMB data analysis. In practice, numerical approximations and the presence of masked 
regions will entail that PI and P2 will only hold approximately; nevertheless, simulations 
have suggested that these properties do ensure a remarkable performance of needlets when 
applied to actual data from CMB experiments, see 



29 



30 



31 



36[ for further developments. 



III. SPIN NEEDLETS 



Throughout this paper, we shall assume that there exist a grid of cubature points {^jk} , 
and a set of corresponding weights {Xjk} such that following discrete approximations of 
spherical integrals hold: 



Ajfc {2^m(0jfc)} {2^£'m'(^ifc)} 



S2 



{2Yem{j)} {2>f'm'(7)}c?7 



(9) 



Here ±2^£m are the spin spherical harmonics defined in the appendix. For the standard scalar 
case, the existence of such points is well-known and provided by many different constructions, 



see for instance 



18 



19| ) and ( |37| ) , see also [20|, l21j for further discussion and 



29 



30 



for extensions to the generalized needlets case. For the spin spherical harmonics we consider 
here, the validity of ([9]) is going to be investigated mathematically elsewhere. We stress, 
however, that even if equation ([9]) is not known, or the points {C,jk} are not explicitly given, 
then one can use any sensible collection of points and weights, and the results will still 
hold approximately; in other words, our results below will continue to hold with minor 
numerical approximations when implemented on any package with a pixelization scheme 
such as HEALPix (this issue is discussed rigorously and in greater detail in [341]). 
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Spin needlets are defined as (see ( 3J]) for a complete matliematical treatment and more 
rigorous results) 

fe(7) = V^E^(:b7) E ^^^^rnm {2^,^(0^)} . (10) 

I m=-l 

A comparison between ([T]) and ( ITOl) highlights immediately that spin needlets make up a 
natural extension of the ideas underlying the approach in the scalar case to the framework 
of a spin field. This deep link between the two constructions should not hide, however, 
some profound differences between ipjk{^) and ipjk;2{'l) as mathematical objects. Indeed, as 
recalled in the previous Section ipjkil) is a standard scalar function which induces a linear 
map ([3]) leading from T(7) —>■ jSjk, i.e. from a scalar quantity to a scalar quantity. On 
the contrary, tpjk;2 induces a linear map leading from spin 2 quantities to spin 2 wavelet 
coefficients. The quantities in ( ITOi) depend on the choice of the coordinate system for 7 and 
C,jk] these two coordinate systems may be chosen independently. If the coordinate system 
for 7 is rotated, ipjk;2{l) transforms like a spin 2 vector at 7, while if the coordinate system 
at is rotated, ipjk;2{'l) transforms like a spin —2 vector at C,jk- '4^jk;2{l) has a precise 
mathematical status (see ([3^])) as a linear map from spin 2 vectors at ^jk to spin two 
vectors at 7. Indeed, if f is a spin 2 vector at ^jfc, ipjk;2{,l)v makes sense as a spin 2 vector 
at 7, since the product of a spin —2 vector and a spin 2 vector at a point is a well-defined 
complex number, independent of choice of coordinates. (Thus it would be more proper to 



write {2Yim{l)} ® \ 2yha{ijk)j than {2Yim{^)} \2yim{ijk)j in (HO]).) 

As we shall show below the fact that ipjk-2{l) is not a scalar function does not prevent use- 
ful applications for the reconstruction and testing on physically meaningful scalar quantities 
such as the angular power spectra Cf^, C^^. We note that in the mathematical results of 



3J], the factor b{l/B^) is replaced by 6(a/ (/ — 2){l + 3)/B^); this reflects the fact that the se- 
quence {(Z — 2)(/ + 3)}j^3 4 represents the eigenvalues of the Laplacian operator associated 
to spin spherical harmonics. However, since our main interest here is high-frequency asymp- 
totics, and since / ~ a/(/ — 2)(/ + 3) for large Z, we use the simpler b{l/B^) here to highlight 
the similarity with the usual presentation of scalar needlets. (Antithetically, one may mod- 



ify the definition of the latter by replacing b{l/B^) with b{\Jl{l + 1)/ B^); note indeed that 
+ 1)}^=! 2 provides the sequence of eigenvalues for the usual Laplacian operator on the 



sphere. We refer to 



29 



30 



3l| for more discussion on this point.) 



In figure [T] we show the projection of '4'jk;2{^) on the plane for B = 1.2 and j = 10. 
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FIG. 1: Projections of ipjk-2{l) for B = 1.2 and j = 10. Left plot: real part, right plot: imaginary 
part, lower plot: modulus. 

We show the real and imaginary component as well as the modulus. In figure [2] we show 
how ipjk;2{l) falls off with the distance from the center for B = 1.2 and j = 10, 20, 30. As 
expected, we clearly see how 4'jk;2i'l) falls of faster for higher values of j and thus picking 
up smaller scales. 

We shall now investigate the localization properties of (|TOl) . Localization properties in 
multipole space are a straightforward consequence of the properties of the weight function 
b{x)] localization properties in pixel space are much less straightforward and require rather 



sophisticated mathematical arguments. The following inequality is established in 3J] and 
generalizes ([7]): 

Proposition 1 [34] For any M = 1,2, ... there exists a positive constant cm such that 
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FIG. 2: The modulus of ipjk-2{l) for B = 1.2 and j = 10 (left plot), j = 20 (right plot) and j = 30 
(lower plot). The angle 6 is the distance from the center point. 



for any point ^ E S"^ we have 



\i^jk;2{l)\ < 



cmB^ 



(l + 5^-arccos(|7-e,fe|))*' 



(11) 



Note that although {iljjk-2{,l)} are not scalar valued, \ipjk;2{,'l)\ is clearly a well-defined 
scalar function, so that the inequalities (fTTl) are consistent. 
Now consider the spin 2 fields 

Q(7) + iU{^) = a'".;2 {2Ylm{l)} 

Im 

where we have introduced the complex-valued random coefficients 



aim;2 = I {Qil) + iU (7)} {2Ylm{l)}d^ = -{aim:E + mm;B) 
1 32 



(12) 
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where aim;E, 0'im;B denote, respectively, the spherical harmonics coefficients of the E,B com- 
ponents of the polarization random field. The spin needlet coefficients are defined as 



: = / {Q{l)+^U{^)}^jk,2{l)d'y 

= a/A^^ &(7r-)a«m;2 {2^m(^iA:)} 



Im 

in obvious analogy to (jl]) for the scalar case, only replacing spin spherical harmonics and 
the corresponding (scalar) random coefficients. Note that 

i m=-e 

is a spin 2 vector at C,jk and a spin —2 vector at 7; (3jk;2 is a spin 2 vector at C,jk- In other 



words, Q{'y)+iU{^) = Y.im^im;2 {2yim{.l)} is spin 2; when muhiplied with V'jA:;2(7), the spin 



2 factor annihilates with the spin —2 factor in 2^m(7); and we are just left with a spin 2 
factor in (3jk-2, given by {2Yim{,ijk)} ■ In figure [3] and H] we show maps of the amplitude and 
direction of polarization of the input map as well as for needlet coefficients at j = 10, 20, 30 
for B = 1.2. 

Assuming the following reconstruction formulae can be shown to hold (see (|34l|) 



{Q(7) +^f/(7)} = $^fete(7) . (13) 

jk 

Note that the right side of f|T3l) makes sense, independent of coordinate system chosen for 
the C,jk, and defines a spin 2 vector at 7. To see this, one need only note again that the 
product of two quantities, one which transforms like a spin 2 vector at C,jk, and the other 
which transforms like a spin —2 vector at C,jk, is a well-defined complex number, independent 
of choice of coordinate system. An example of the reconstruction is shown in figure [5l In 
the figure we show the relative difference in percent between the input and reconstructed 
Q and U maps. Outside the masked regions, the difference is equal to the numerical noise 
obtained by doing a simple spherical harmonic transform followed by the inverse transform. 
Close to the masked regions, there is a small area where the reconstruction error is larger. 
The size of this region will depend on the choice of B. In this case B = 1.2 which is a rather 
low value resulting in a relatively large region with larger error around the mask. Note that 
most of these points still have reconstruction errors smaller than 5%. We will discuss the 
accuracy of reconstruction in detail in a following paper. 
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input map, amplitLide 6',, arnpliiijde, | = 10 




^jt, amplitude, j=20 




FIG. 3: The input CMB map as well as the needlets coefficients for B=1.2 and j=10, 20 and 30. 
The plots show the polarization amplitude. 



As detailed in the previous Section, a capital property for the random needlet coefficients 
in the scalar case is asymptotic uncorrelation at any fixed angular distance (for a smooth 
angular power spectrum), as the frequency parameter diverges ( j oo). A natural question 
is the extent to which this property may continue to hold in the spin case. The answer to 



this question is provided in [341 , where it is shown that, under mild regularity conditions on 



the angular power spectra, for all M > there exist positive constants Cm such that 

I {Pjk;2f3jk';2) \ ^ Cm_ , . 

(IfeO - (l + i?^-arccos(|Ofc-0.'|))*^ ■ ^ ^ 

Again, \{Pjk;2Pjk';2)\ is a well defined scalar functions, despite the fact that {l3jk;2l3jk';2) 
are spin quantities that depend on the choice of coordinates for tangent planes. (|T^ sug- 
gests that the spin spherical needlet coefficients can be consistently used for angular power 
spectrum estimation and map reconstruction, as detailed in the following section; in the 



scalar case, this argument is rigorously derived in ( 20|]), whereas for the spin situation we 



are concerned with here, we refer again to 3^] for a more complete mathematical analysis. 
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FIG. 4: The input CMB map as well as the needlets coefficients for B=1.2 and j=10, 20 and 30. 
The plots show the polarization direction. 
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FIG. 5: The relative difference (in % between an input CMB map (Q and U map) and tlie 
reconstructed map in the presence of a mask (in this case the P06 galactic cut used by the WIVlAP 
team for polarisation). 
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Now assume that we have a masked region M where no polarization data are actually 
available; of course, this implies that we shall only be able to recover the coefficients 



{Q{^)+ium{2yimmdi ■ 



(15) 
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Due to the poor localization properties of spherical harmonics, the coefficients {aim;2} need 
not be close in any meaningful sense to {aim;2} ', consequently any statistical procedure based 
naively upon them (including estimation of Cf"^, or the reconstruction of the scalar 
maps -^(7), -6(7)) is likely to be seriously flawed. On the other hand, consider 



I 



Im 



for locations {^jk} that are "sufficiently away" from the masked region M, i.e. d{^jk, M) > 
5 > 0, it is easy to see that (see ffTTl) ) 



Pjk;2 ~ f3jk;2 



{g(7)+^f/(7)}te(7)c^7 



M 



< 



cmB^ 



(l + 5Jarccos((5))^'^ 



\{Q{i)+tum\di 



M 



(16) 



i.e. the coefficients Pjk;2, Pjk;2 become asymptotically equivalent at high frequencies j 00. 
It should be noted that throughout this paper we have decided to adopt the spin spherical 



harmonics formalism considered by 



39| . A completely analogous result could have been 



obtained by taking combinations of b{l/B^) with grad and curl harmonics, i.e. starting 



from the approach of 



401 ]. We do not focus on this possibility here for brevity's sake. 



IV. SOME STATISTICAL APPLICATIONS 



The purpose of this Section is to provide some examples of the statistical applications 
which can be entertained by means of spin needlet coefficients. Our purpose here is not 
to provide recipes which are ready-to-use for CMB data analysis - this certainly requires 
much more computational and analytical work to take into account the features of CMB 
maps: foregrounds, anisotropic noise, multichannel observations and many others. Our 
purpose here is different: it is well-known how wavelets in the scalar case have proved to be 
a valuable instrument when dealing with a variety of data analysis issues. The discussion 
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below suggests that spin needlets can be just as important. Due to their specific nature 
as mathematical (spin) objects, these applications however require extra care, even in the 
idealistic circumstances we focus on here. 



A. Estimation of the sum of the angular power spectra Cf^ + Cf^ 

We recall that an estimator of the binned temperature power spectrum can be constructed 
from scalar needlets, as follows. Take 



def 



in the idealistic case of no mask, it is readily seen that 

/ 



T, 



J2 b'i—){2l + l)Q 

B3-^<l<Bi+^ 



. r, 



In the idealistic circumstances where the sky is fully-observed, the previous sum includes 
Nj terms, where we recall that Nj is the number of pixels at the resolution j. The result 
continues to hold (with some corrections to take into account the fraction of sky coverage) in 
the presence of a sky cut, although clearly Ylk runs over a smaller number of terms (in other 
words, a sky- fraction correction factor must be introduced). This estimator was proposed 



m 
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22l|): in 221] applications are also shown to the cross-correlation between CMB and 



Large Scale Structure data, in 27|, 



of observational noise, whereas in 
and asymmetries in CMB maps. In 



this estimator is extended to allow for the presence 
the same approach is applied to search for features 



201 1) it is also shown that, in the Gaussian case 
-diV(0,l), 



Var { Tj 



(17) 



where — >d denotes convergence in distribution and A^(0, 1) a standard Gaussian law. A result 



like (fT7|) provides also a clue for testing goodness-of- 
high frequencies) for the angular power spectra. In 



it and driving confidence intervals (at 



2l| the behaviour of statistics such as 



Tj are considered for partial regions of the sky, also as a device for testing asymmetries. 
From the mathematical point of view, results like ( fT7|) are entirely justifiable on the basis 
of the uncorrelation properties we discussed earlier in P1-P2 (see ([H])). As the same sort of 
uncorrelation property has been established in (fT^ . it is natural to investigate whether a 
similar procedure for the estimation of the binned angular power spectra is feasible here. 
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The answer turns out to be positive, in the following sense. Consider the (scalar valued) 
estimator 



i;2 



Starting from the idealistic case of no sky-cuts, we obtain easily 



Im 



def 



j;2 , 



because 



m 

and as before we took A^^ = Atc/Nj. Here 



2/ + 1 
An 



Ti —< \aim-2f > — < \aimE + idimsf >— Cf^ + Cf^. 



Likewise, it is possible to show that 



m 



rj;2 - rj; 2 
Var \ r,;2 



>d A^(0, 1) , as j ^ oo , 



(18) 



i.e. it is possible to prove that is consistent and asymptotically Gaussian around its 
expected value Tj.2. Note that, as the electric and magnetic components of the polarization 
field are uncorrelated, the variance in the denominator is simply the sum of the variances of 



2G 



2l|. 



the two scalar components. Likewise, the asymptotic theory can be developed as in 
As mentioned before, the presence of masked regions of the sky requires the introduction 
of a sky coverage fraction. The presence of anisotropic noise is more interesting and can 



be dealt with along the lines of 



281]), i.e. by introducing weighted, rather than simple. 



averages of the squared needlet coefficients. 

The previous procedure allows one to estimate (a binned form of) the angular power 
spectra Cf^ + C^^^. Although this could be sufficient for some purposes, it is clear that 
for CMB applications we are interested, in general, in the estimation of Cf^ and C^^ as 
separate quantities; this is an issue to which we shall come back later in this Section. 
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B. Testing for non-Gaussianity 



As a further statistical application, it is possible to consider the investigation of non- 
Gaussianity in the joint law of temperature and polarization data. A standard idea to 
focus on non-Gaussianity is to consider the skewness and kurtosis of wavelet coefficients; 
for brevity, we shall concentrate on the latter statistic. To normalize our coefficients, we 
estimate their variance by 



where we define also 



= < > = 



> 



As a consequence of (fT8|) and along the same lines as in 20||, it can be shown that 



p lim 



a 



p lim — 



plimj^oo denoting as before convergence in probability; loosely speaking, this is to say that 
the mean of the squared spin needlet coefficients is a consistent estimator of their variance. 
We focus then on the normalized coefficients 



'jk;2 



def (3jk;2 



(Jo 



to test for the joint Gaussianity of the temperature and polarization fields, we may consider 
for instance their kurtosis, i.e. 



Kj.pp 



E |3 



Jjk;2 



which converges to zero in probability as j ^ oo. More interestingly, we might be interested 
in testing for the joint Gaussianity of polarization and temperature maps. Likewise, it seems 



possible to focus on the needlet 
suggested for the scalar case by 



jispectrum for joint temperature and polarization data, as 
25|. These issues will be developed in a future paper. 



C. Reconstruction of the E and B maps 

A common strategy to reconstruct the (scalar) maps of the electric and magnetic (or 
grad and curl) components of the polarization field is well-known; starting from the maps 
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{Q ± ill} , the scalar coefficients are evaluated by means of the Fourier inversions (fT2l) . The 
maps are then recovered by the standard spectral expansion. As simple as this procedure 
can seem, it is prone to severe problems in the analysis of actual data, where the presence of 
masked regions can lead to severe errors in the evaluation of exact Fourier transforms. We 
shall show here how one can use needlets to produce scalar maps by a novel technique. 

We recall ffist how to obtain scalar maps from polarization data; indeed, as in we 
focus on 



E(7) 



(Q+tum + {df{Q~tum 

where we have used the spin raising and spin lowering operators {d, d) defined in the ap- 
pendix. The crucial property to recall is 



{i-2) 
'(^ + 2) 



{i-2)\ 



Yemil) 



Hence we obtain, in view of (fT3|) 



E(7) 



TT\ 2 



jk jk 



jk 



i / (f + 2)! 
B^'V {£-2)1 



jk 



im 



B^^y {£-2)1 



{Yem{l)} {2YU^jk)} 



{Yem{l)} {2YeU^jk)} 



{f^jk;2Vjk{l) + l3jk-2Vjk{l)} , 



jk 



where 



im 



TT\2 



(19) 



(20) 
(21) 
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which is scalar in 7, spin —2 in ^j^. Likewise 



5(7) 



(Q + ^Um-{^)'{Q-^Um 



jk jk 



^ E/^^'^;^ J-) {Yemm {2YeU^,k)} 



jk 



im 



jk 



£ , /(£ + 2)! 



im 



{Yim{l)] {2Y,m{ijk)} 



Y {(^3k;2^jkin) - I3jk,2^jk{l)] ■ 



jk 



In the presence of fully observed maps, it is immediate that 

jk K Im ) im' 



I ^ l ii + 2)\ 



(22) 



{Yim'm {2Yemi^jk)} 



^ E ^ (E ^kb{^)-^m;2 {2YUi3k)) ] ^ 

jk K Im J im' ]j \ )■ 



{Yem'il)} {2Yem'{^jk)} 



I E("'-'2 + a/-m;2)>"^m'(7)\/ 1^3^ E ^^^^^^^^^^^^ 



Im 



ml 

m 



} , aim;E \ I —„ yim 17 J 



Im 



ii-2) 



(23) 



likewise, it is immediate that 



^(7) = E^''"'-si 



Im 



(1 + 2)1 

{l-2)\ 



YUl) , 



(24) 



as expected. This is just a rephrasing of f|T3l) in terms of the underlying [electric and 
magnetic) scalar fields. 

In the presence of masked maps, E and B are unfeasible; our idea is then to use thresh- 
olding techniques, which have been shown to be powerful when combined with the wavelet 
approach (see for instance 43], or 2^ for applications in a CMB framework). In particular, 
assume a portion M of S'^ is masked; for each pixel {^jk} we can define the fraction 

/s2/M lte(7)l^^7 



Sjk 



Js2 \ll^jk;2{l)f 
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This ratio is clearly measuring the amount by which the needlet coefficient localized at {^jk} 
is corrupted by the presence of missing observations. Note that, due to the localization 
properties of spin needlets, Sjk will converge to unity for all pixels {^jk} which are outside 
the masked regions. We can now define the thresholded parameters 

Pjk;2 = Pjk;2^{Sjk > tj) , 

where the function I{sjk > tj) takes the value one if Sjk > tj, zero otherwise; tj is a 
thresholding parameter which is assumed to converge to one as j oo. We can then 
consider the reconstructed maps 

^*(^) + ^^(7)] I , (25) 

^lY. fe;2^J-fc(7) - :^^(7)] • (26) 

jk 

In view of f|T6l) . we expect E*{'y) ~ -^(7), B*{'j) ~ -8(7); the validity of these approximations 
depends upon the distance of 7 from the masked region and the cosmic variance (i.e., the 
lower the value of jCf^^, C;^^} at low multipoles, the better the approximation). The 
accuracy of this approach will be tested extensively in a future publication. 

It should be noted that for the E and B fields in ( l23l we used the same definition 
and notation as given by ([42]); this definition for the scalar components seems somewhat 
natural as it follows directly from the application of the spin raising-spin lowering operation 
to polarization data. More commonly, in the CMB literature the E and B random fields are 
defined instead as 

^(7) = J^a«m;£;>^m(7) , ^(7) = a/m;B>^«m(7) ; 

Im Im 

i.e., these fields differ by a normahzation factor a/(/ + 2)!/(/ — 2)! multiplying the random 
spherical harmonic coefficients {ai^.E, «/m;B}- It is indeed possible to obtain needlet approx- 
imations of these maps, by focussing on 

I jk 

jk 
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where 

^Ai) = v^E^(:^) ^^^-(^)> i^fe-^)} • (27) 

im 

It should be noted however that E{^), -8(7) do not follow trivially by application of the spin 
raising - spin lowering operators to polarization data; likewise, and as opposed to (!20!) . (127|) 
is not the outcome of differentiation by d (d) on spin needlets, so that its properties are less 
clear. (Because of (pTj) . (1201) is well-localized by the results of |34|.) 

The results in the present subsection suggest that spin needlets may allow for new sta- 
tistical techniques when dealing with map reconstruction, even besides those allowed by 
wavelet methods in the scalar case. More generally, due to the poor localization properties 
of standard spherical harmonics, in the presence of masked regions and instrumental noise 
the construction of E and B from the coefficients {0,1^-2} (see (ITSl) ) may not ensure satisfac- 
tory results (recall {aim;2} need not be close to {aim;2} in any meaningful sense for partially 
observed polarization maps). On the other hand, in view of (JT6l) the coefficients |/3jfc| , are 
much less affected by masked regions, despite the fact that these quantities can themselves 
be viewed as linear combinations of the spherical harmonic coefficients. This suggests that 
spin needlet coefficients may be extremely useful when attempting to build optimal polar- 
ization maps; for instance, it seems natural to suggest the implementation of Internal Linear 
Combination techniques based on wavelet coefficients (as done in the scalar case by 26|). 
without the need to go through scalar maps as a first step. These developments are left for 
future work. 



D. Estimation of the angular power spectra Cf^,Cj^^ 

We are now ready to address the estimation of the angular power spectra Cf^, C^^ as 
separate quantities, which is clearly much more meaningful from a physical point of view. 
For this purpose, it suffices to write 



Eji'f) = -lj2^'i^)iaim;2 + ai-m;2)YlUl) 



Im V \ ' 



2^b {^)al„^■,EYMJ———^^ , 

Im V V / 
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and likewise 



Clearly 



Im V ^ ' 



< / ^;(7)ci7 >= / < ^- (7) > dn 



< 



^BB 
' j 



Et2( ^ \nEE + 2)! def ~ee 

In the presence of a mask, it is then natural to suggest the following (unbiased) estimator: 

here, the denominator is just a sky-fraction normalizing factor, whereas the integrals can be 
easily approximated by sums in pixel spaces. It is obvious that these estimators are unbiased 
for r^^,r^^; likewise, one might focus on 

which are unbiased for the binned spectra 'Ylii^'^i^ / B-')Cf^ ^'Y^i}P'{l / B^)CP^ . A more com- 
plete set of properties and implementation on polarization data will be provided elsewhere. 



V. A COMPARISON WITH ALTERNATIVE CONSTRUCTIONS 

To compare our approach with possible alternative constructions (see also f^), let us 
consider again the scalar fields of E- and S-modes, which we recall can be written as 

-^(7) := ^^Cilm;EYlm{l) , ^(7) := ^ a;^;^!^/™ (7) . 
Im Im 

In the idealistic case of fully observed and noiseless maps for -^(7) = {Q{l) +*f^(7)} , the 
coefficients {aim-^E^O'im;B} could be exactly recovered and it would be natural to derive the 
scalar random needlet coefficients on E{j), B{j) as 
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In other words, one might propose to define a wavelet construction by deriving first the 
scalar maps, and then applying on them the standard procedures. The purpose of this 
Section is to compare this approach with the one we suggested earlier to show that they 
are not equivalent. More precisely, we shall address the question about the localization 
properties of (125]) . and investigate the effect of the presence of a masked region M. 

Recall that for scalar fields, when is "far enough" from M we have the approximation 



Im ^ 



In the presence of a masked region we obtain for instance 



l^jk - ^/^^b (—j aim-EYlmi^jk) 
Im, ^ 



2 

Now note that 



afm = -U I [^(^) (2^'- (7)) + ^(7) {-2Ylm (7))] ^7 

S2/M 



^ - Pfk = -V^kY^b ^ {at - afj Yim {ijk) 

Im ^ 

-^^E^ (iy) {/^ [^(7) (^(7)) +m (^(7))] d^Y'^ii^'^ 

[ni)X,k;2 (7) + ni)X,k;-2 {l)] dj , (29) 



where 



Xjk;2 (7) = V^E ^ (2^™ (7)) Ylm (Ofc) , (30) 

Xjk;-2il) = ^/>^kY,b(^^UYlrnmy^i^jk) ■ (31) 

Now the properties of fl30l) . fl3T|) are still unknown and lack a rigorous mathematical 
exploration. Indeed the function J2m (^2^Zm (7) j Yim {^jk) does not represent any form 
of projector on a proper subspace, as it is the case for the corresponding expression 
J2m (^im (7) j {^jk) fof the Standard needlet construction in the scalar case. In fact 
looking at figure [6] (compare with figure [2] for the spin needlets) we see that the shape of 
\Xjk;2 (7) I is such that at the origin the function has a zero, rather than a maximum, and 
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FIG. 6: The modulus of Xjk;2{'l) for B = 1.2 and j = 10 (left plot), j = 20 (right plot) and j = 30 
(lower plot). The angle 6 is the distance from the center point. 



then increases away from the center. In other words, while spin needlets enjoy a number of 
properties which are analogous to those in the scalar case (good approximation properties 
in the presence of a mask (see ( JT6l) above), tight frame property, asymptotic uncorrelation), 
these properties are not in general known to be shared by constructing wavelets on the scalar 
random fields E, B after implementing the spherical harmonic transforms, in the presence 
of masked regions. 

We believe these remarks provide a clear rationale for our construction of spin needlets, 
which are directly embedded into the geometric structure of polarization random fields. 
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VI. APPENDIX: SPIN SPHERICAL HARMONICS 



In this Appendix, we review a few basic facts about the analysis of polarization data. Our 



review will be quite informal, and we will follow mainly the formalism by [39|; we refer also 
to Qf„.a„e,.va.e.ap,„ac. andante, more detailed treatment, see also [42|. All 
the material discussed here is absolutely standard, and it is reported just for completeness. 
A standard route for the expansion of polarization random fields is based on spin spherical 



harmonics (see for instance 



39|,[42|)). Consider the standard representation of the group 



of rotations 5*0(3) by means of the Wigner's D-matrices, with elements 



t=OV{m-n) 

where we recall that 
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The spherical harmonics of spin n can then similarly be defined as 



2/ + 1 r^U 



(<^,?9,0) 



'21 + 1 



(32) 
(33) 



An important property we shall repeatedly use is 



see for instance 



m=—l 



21 + 1 
Air 



(34) 



34| . Theorem 2.8. 



Under a rotation around the tangent plane at the point (t?, (p), „Y;m(7?, ip) transforms into 
e~'^^^Yim{'& , ^p) , i.e. as a spin n function. Spin n functions can be defined iteratively by 
means of the spin raising and lowering operators 



We have for instance 
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d'd sin d d(p 



sin"'(? 



(nG(^,^)) 
(nG(^9,^)) 



(35) 
(36) 



d [nYimi'd, ip)) = [{l-n){l + n + UiYimi'd, y^)) 
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whence 

which generahzes iteratively to (for < n < /) 



{i + ny. 



and (for > n > —I) 



^Ylm{^, if) 



{l + n)\ 
{l-n)\ 



1/2 



More exphcitly 



2Yl„i{'d,Lp) 



(/-2)!1'/' 



L(^ + 2)! 

(Z-2)!]'/' 
(7T2)! 



sm v — H — : — - t:— I sm v 



d'd sin dip 



d id 

+ — 



dd sin dip 



— sm 17 I H — ^ — - -TT- I sm v 



d-& sin ?9 



TTl d 

YU^,v)-—Yi^{d,^) 
sm -(7 dv 



Simple computations yield 



(^ + 2)! 
{l-2)\ 



{±2^m} 



m 



sin^ sin 'd 



2m cot -(9 d d"^ 

+ 2m ^ — - — cot "i/TT-r + 



sint? 



d'd d^"^ 



Y 



Im ■ 



(37) 



We can then introduce the E and B scalar components, whose random coefficients coefficients 
can be obtained (see 



39|],j42|)) evaluating the spin 2 transforms 



±2Qlm-= I Q{x) {±2Ylm{x)} dx , ^2Ulm ■= f/(x) {±2^m (a^) } C^a; , 

'52 Js^ 



and then proceeding to evaluate the coefficients 

^2<5/m + hUlrn^ + ^-2<5«m — i-2Ulm 
2Qlm + hUlm ) — ( -2Qlm — i-2Ulm 



(38) 



(39) 
(40) 
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